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By making use of the Lewis-Riesenfeld invariant theory, the solution of the Schrodinger equa- 
tion for the time-dependent linear potential corresponding to the quadratic-form Lewis-Riesenfeld 
invariant I q (t) is obtained in the present paper. It is emphasized that in order to obtain the gen- 
eral solutions of the time-dependent Schrodinger equation, one should first find the complete set of 
Lewis-Riesenfeld invariants. For the present quantum system with a time-dependent linear poten- 
tial, the linear I\(t) and quadratic I q (t) (where the latter I q (t) cannot be written as the squared of 
the former I\(t), i.e., the relation I q (t) = cl 2 (t) does not hold true always) will form a complete set 
. of Lewis-Riesenfeld invariants. It is also shown that the solution obtained by Bekkar et al. more 

' recently is the one corresponding to the linear I\(t), one of the invariants that form the complete set. 

In addition, we discuss some related topics regarding the comment [Phys. Rev. A 68, 016101 (2003)] 
O ■ of Bekkar et al. on Guedes's work [Phys. Rev. A 63, 034102 (2001)] and Guedes's corresponding 

reply [Phys. Rev. A 68, 016102 (2003)]. 
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I. INTRODUCTION 



Recently, Guedes used the Lewis-Riesenfeld invariant formulation [1] and solved the one-dimensional Schrodinger 
equation with a time-dependent linear potential [2]. More recently, Bekkar et al. pointed out that [3] the result 



obtained by Guedes is merely the particular solution (that corresponds to the null eigenvalue of the linear Lewis- 
Riesenfeld invariant) rather than a general one. In the comment [3], Bekkar et al. stated that they correctly used the 
invariant method [1] and gave the general solutions of the time-dependent Schrodinger equation with a time-dependent 
linear potential [3]. However, in the present paper, I will show that although the solutions of Bekkar et al. is more 
general than that of Guedes [2], what they finally achieved in their comment [3] is still not the general solutions, either. 
On the contrary, I think that their result [3] also belongs to the particular one. The reason for this may be as follows: 
according to the Lewis-Riesenfeld invariant method [1], the solutions of the time-dependent Schrodinger equation can 
be constructed in terms of the eigenstates of the Lewis-Riesenfeld (L-R) invariants. It is known that both the squared 
of a L-R invariant (denoted by I(t)) and the product of two L-R invariants arc also the invariants, which agree with 
the Liouville-Von Neumann equation -jLl(t) + \ [I(t), H(t)] — 0, and that if I a and h are the two L-R invariants 
of a certain time-dependent quantum system and \ip(t)) is the solution of the time-dependent Schrodinger equation 
(corresponding to one of the invariants, say, I a ), then If,\ij;(t)) is another solution of this quantum system. So, in 
an attempt to obtain the general solutions of a time-dependent system, one should first analyze the complete set of 
all L-R invariants of the system under consideration. Historically, in order to obtain the complete set of invariants, 
Gao et al. suggested the concept of basic invariants which can generate the complete set of invariants [4] , as stated 
in Ref. [4], the basic invariants can be called invariant generators. As far as Bekkar et a/.'s result [3] is concerned, 
the obtained solutions are the ones corresponding only to the linear invariant (i.e., I\(t) = A(t)p + B(t)q + C(t)) 
that is simply one of the L-R invariants, which form a complete set. It is apparently seen that the quadratic form, 
7 q (t) = D(t)p 2 + E(t)(pq + qp) + F(t)q 2 + A'(t)p + B'(t)q + C'(t), is also the one that can satisfy the Liouville-Von 
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Neumann equation, since it is readily verified that the Lie algebraic generators of I q (t) form a Lie algebra, which 
possesses the following commutators 

[q 2 ,p 2 ] = Mpq + qp), [m + qp, q 2 } = -4i<? 2 , [pq + qp,p 2 ] = &p 2 , 

[q,p 2 ] = 2ip, [p, q 2 ] = -2iq, [q,pq + qp] = 2iq, [p,pq + qp] = -2ip. (1.1) 

However, for the cubic-form invariant, it is easily seen that there exists no such closed Lie algebra. This point holds 
true also for the algebraic generators in any high-order L-R invariants 7™. So, it is concluded that for the driven 
oscillator, only the linear I\(t) and quadratic I q (t) will form a complete set of L-R invariants. Note that here I q (t) 
should not be the squared of I\(t), i.e., I q (t) ^ cl 2 (t), where c is an arbitrary c- number. It is emphasized here that 
Bckkar et aVs solution is the one constructed only in terms of the eigenstates of the linear invariant I\(t). Even though 
only for the linear invariant I\(t) Bckkar et al.'s result [3] can truly be viewed as the complete set of solutions, it still 
cannot be considered general one of the Schrodinger equation, since the latter should contain those corresponding 
to the quadratic invariant Iq(t). In brief, Bekkar et al.'s solution and my solution, which will be found in what 
follows, together constitute the complete set of solutions of the Schrodinger equation involving a time-dependent 
linear potential. 



II. ON THE COMPLETE SET OF L-R INVARIANTS 



According to the L-R invariant theory [1], if the eigenstate of the linear invariant Ii(t) corresponding to A„, i.e., 
one of the eigenvalues of h{t), is |A„, t), then the solution of the Schrodinger equation can be written in the form 



|*(£))schr = J^c„ exp 



-Mt) 



\\n,t) 



(2.1) 



with c n 's and 0„(i)'s being the time-independent coefficients and time-dependent phases [1,4], respectively. This, 
therefore, means that the solutions of the time-dependent Schrodinger equation can be constructed in terms of the 
complete set of eigenvector basis, {|A„,t)} ; of I\(t). Moreover, one can readily verify that the squared, I 2 , of the 
linear invariant is the one satisfying the Liouville-Von Neumann equation, and that /i(i)|^E'(i))schr is also a solution 
(but not another new general one) of the same time-dependent Schrodinger equation, since it is readily verified that 
/i(t)|^ / (i))schr can also be the linear combination of the eigenstate basis set {|A„,t)} of I\{t), i.e., 



/i(i)|*(i))schr = ^6„exp 



I Art) t)i 



(2.2) 



where the time-independent coefficients 6„'s are taken b n — X n c n , which is obtained via the comparison of the 
expression (2.2) with (2.1). 

Thus, the above discussion shows that the linear invariant I\ and its squared I 2 have the same eigenstate basis set 
and therefore I\ and I 2 cannot form a complete set of L-R invariants. In contrast, if for any c-number c, the quadratic 
I q cannot be written as the squared of linear I\ with various integral constants A , B and Co (for the definition of 
Aq, Bq and Co, see, for example, in Ref. [3]), namely, the relation I q = cl 2 is always not true, then {Ii,I q } is the 
complete set of L-R invariants, which enables us to obtain the general solutions (complete set of solutions) of the 
time-dependent Schrodinger equation. 

Perhaps someone will ask such question as, "Does there really exist such quadratic I q that can always not be written 
in the form d 2 T or "Maybe any I q that satisfies the Liouville-Von Neumann equation can surely be written as the 
squared of certain I\. Really?" Now I will discuss these questions. Consider a given quadratic invariant I q that 
is written I q (t) = D(t)p 2 + E(t)(pq + qp) + F(t)q 2 + A'(t)p + B'(t)q + C'(t) whose time-dependent parameters are 
determined by the Liouville-Von Neumann equation, and a certain linear invariant I\(t) = A(t)p + B(t)q + C(t), the 
squared of which is I 2 = A 2 p 2 + AB(pq + qp) + B 2 q 2 + 2C(Ap + Bq+ §). Since the functions A, B and C can also 
be determined by the Liouville-Von Neumann equation, the only retained parts left to us to determine is the integral 
constants A , B and C . Choose the appropriate integral constants in A, B and C, and let I q be the squared of I\ 
(should such case exist), and then we have 



D = cA 2 , E = cAB, F = cB 2 , 
A' = 2cAC, B' = 2cBC, C = cC 2 . 



(2.3) 
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If a given / q can really be written as the squared of I\, the above six equations are just used to determine the c- 
number c and the suitable integral constants A n , B and Co in the functions A, B and C. It is seen that there are 
only four numbers expected to be determined, and that, in contrast, we have six equations. So, it is possible that 
there exist potential parameters c and Aq, Bq, Co which will not agree with Eqs.(2.3) always for a given parameter 
set {D,E,F,A',B',C}, or, for a given parameter set {D, E, F, A' , B' ,C'} there are always no such parameters c 
and Ao, B , Co which satisfy Eqs.(2.3). The existence of I q that cannot be written as the squared of any I\ is thus 
demonstrated. 

So, in the above we indicate that such two invariants I\ and I q (which are independent) form a complete set of L-R 
invariants. 



III. UNITARY TRANSFORMATION ASSOCIATED WITH L-R INVARIANTS 

Now I will solve the time-dependent Schrodinger equation, of which the time-dependent Hamiltonian [2] is given 

H(t) = ^ + f(t)q, (3.1) 

by making use of the Lewis-Riesenfeld invariant theory [1]. The time-dependent L-R invariant used here takes the 
form 

J q (t) - D(t)p 2 + E(t)(pq + qp) + F(t)q 2 + A(t)p + B(t)q + C(t). (3.2) 
With the help of the Liouville-Von Neumann equation, one can arrive at 

2E F 
D+ — =0, E+ — = 0, F = 0, 

m m 

A+ 2Df = 0, B-2Ef = 0, C-fA = (3.3) 

m 

with dot denoting the derivative with respect to time t. The above six equations (referred to as the auxiliary equations 
[4]) can be used to determine all the time-dependent parameters A(t), B(t), C(t) and D(t), E(t), F(t). 

In accordance with the L-R theory, solving the cigcnstates of the invariant (3.2) will enable us to obtain the 
solutions of the time-dependent Schrodinger equation. But, unfortunately, it is not easy for us to immediately solve 
the eigenvalue equation of the time-dependent invariant (3.2), for the invariant (3.2) involves the time-dependent 
parameters. So, in the following we will use the invariant-related unitary transformation formulation [4], under which 
the time- dependent invariant in (3.2) can be transformed into a time-independent one ly, and if the eigenstates of ly 
can be obtained conveniently, the eigenstates of I q (t) can then be easily achieved. 

Here we will employ two time-dependent unitary transformation operators 

Vi(t) = expfo(t)g + 0(t)p], V 2 (t) = cxp[ a (t)p 2 + p(t)q 2 } (3.4) 

to get a time-independent ly . The time-dependent parameters r\, (3, a and p in (3.4) are purely imaginary functions, 
which will be determined in the following subsections. Since the canonical variables (operators) q and p form a non- 
semisimple Lie algebra, here the first step is to transform I q (t) into h(t), i.e., Ii(t) = V^{t)I q {t)V\{t), which no longer 
involves the canonical variables q and p, and the retained Lie algebraic generators in I\{t) are only p 2 , pq + qp, q 2 . 
Note that these three generators also form a Lie algebra (see the commutators (1.1)) . The second step is to obtain 
the time-independent ly, which will be gained via the calculation of ly = V% (t)Ii(t)V2(t). In this step, the obtained 
ly has no other generators (and time-dependent c- numbers) than p 2 and q 2 , namely, ly may be written in the form 
ly = q{p 2 + q 2 ) with q being a certain parameter independent of time. It is well known that the eigenvalue equation 
of ly is of the form Iy\n, q) = (2n + l)<r|n, q), where \n, q) stands for the familiar harmonic-oscillator wavefunction. 
Hence, the eigenstates of the time-dependent L-R invariant I q (t) in (3.2) can be achieved and the final result is 
V\{t)V2{t)\n, q) with the eigenvalue being (2n + !)<;. 



A. The calculation of h(t) = (t)I^(t)Vi(t) 



By the aid of the Glauber formula, one can arrive at 
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h(t) = Dp 2 + E(pq + qp) + Fq 2 + [A + 2i(E/3 - Dn)]p + [B + 2i(F0 - E v )]q 

+ C- [-i(B(3 - Arj) + Dt] 2 + F(3 2 - 2Ej3rj}. (3.5) 

If the two relations 

A + 2i(Ef3 - Drj) = 0, B + 2i(Ff3 — Er}) = (3.6) 
are satisfied, then we can obtain 1 

h(t) = D(t)p 2 + E(t)(pq + qp) + F(t)q 2 . (3.7) 
It follows from (3.6) that the time-dependent parameters in the unitary transformation V\(t) are expressed by 

EB-FA DB-EA 

V ~ 2i(E 2 - DF) ' P~ 2i(E 2 -DF)' ( ' 



B. The calculation of I v = V%(t)h(t)V 2 (t) 

By using the Glauber formula, one can arrive at 

I v = V}(t)h(t)V 2 (t) = Vp 2 + £(pq + qp) + Tq 2 , (3.9) 
where V, £ and T are of the form 

_ ^ AiEa , , N i — 8(Fa — Dp)a r , ,„„ N i 
V = D + r sinh 16pa 2 + ^ ^l— cosh(16pa 2 _ 1 

2i(Fa — Dp) . , .„ n .1 _ , .1 

£ = — - r ^sinh(16pa)5 + E cosh(16 pa) ? , 

(16 pa) 2 

T = F + sm h(16pc^ + 8 ^ Fa D p)p [cosh(16p a ) - ll , (3.10) 

(16pa)i 16pa L J 

respectively. It follows that if the following two equations are satisfied, 

S = Csinh(16/9a)5, 2l ( Fa -^P) = _^ cos h(16pQi)^ (3.11) 
(16 pa) 2 

then the coefficients of pq + qp in Iy is vanishing. In order that we can analyze the above equations (3.11) conveniently, 
the time-dependent parameters a, p (which are expected to be determined) and F, D are respectively parameterized 
to be 

u9 v9 

a = — , p = — , F = h cosh(^/mj9) , D = g cosh(y / uu0). (3-12) 
Substitution of the expressions (3.12) into (3.11) yields 

£ = CsinhU»), i{k "Zl v) =-C, (3.13) 
IJuv 



which can determine £ and 9 (expressed in terms of E, h, g and u, v). It is noted that if the functions u and v are 
finally determined, then the time-dependent parameters a and p in the unitary transformation operator V^(f) (3.4) 
can be obtained. 



1 In general, for the case of three-generator Hamiltonian (the generators of which form a non-semisimple algebra), the time- 
dependent c-number C(t) — [—i(B(3 — Arf) + Dr/ 2 + F/3 2 — 2EPrj\ in h(t) are vanishing. See, for example, in Ref. [5], which is 
a special case of the present problem. 
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In what follows wc will determine u and v via setting T> — T = <r with <; being constant (i.e., time- independent). 
Insertion of (3.12) into (3.10) will yield 

D + kU ~ g ^ [cosh(v^g) - 1] = g, F - ^ ~ g " [cosh( V^fl) - 1] = (3.14) 

Eq.(3.14) can determine the functions u and u, although the problem is very complicated. Here it should be noted 
that 9 which has been determined by (3.13) is also the function of u and v. Thus, in principle, we can obtain the 
time-dependent functions a and p in the second unitary transformation operator V 2 (t) = exp[a(t)p 2 + p(t)q 2 ]. 

Now under the unitary transformation V\(t)V 2 (t) the time-dependent invariant J q (t) is changed into a time- 
independent one, i.e., 

I V = WMViit)]* k(t) [Vi(t)V 2 (t)} = ,(p 2 + q 2 ) (3.15) 

whose eigenvalue is (2n + 1)<T and the corresponding eigenstate is |n, q) that is the familiar stationary harmonic- 
oscillator wavefunction, and the eigenvalue equation of the time-dependent invariant I q (t) is thus given as follows 

I q (t)Vi(t)V 2 (t)\n,q) = (2n+l)<;V 1 {t)V 2 (t)\n,q). (3.16) 



C. The solutions of the time-dependent Schrodinger equation 

According to the L-R invariant theory, the particular solution \n, t} Schr of the time-dependent Schrodinger equation 
is different from the eigenfunction of the invariant I q (t) only by a phase factor exp [i</> n (i)] , the time-dependent phase 
of which is written as 

4> n {t) = f\n,q\ [V 1 (t')V 2 (t')}' ! [H(t')-id/dt'] [V^V^t')] \n,q)dt'. (3.17) 
Jo 

This phase <j> n {t) can be calculated with the help of the Glauber formula and the Baker-Campbell-Hausdorff formula 
[6,7]. 

The particular solution \n, i) gchr of the time-dependent Schrodinger equation corresponding to the invariant eigen- 
value (2n + 1)<; is thus of the form 

V^t)V 2 {t)\n,q). (3.18) 
Hence the general solution of the Schrodinger equation can be written in the form 

l*(?,i)}schr = $>„|M> Schr , (3.19) 

n 

where the time-independent c-number c„'s are determined by the initial conditions, i.e., c n =schr (n,t = 0\^(q,t = 
0))schr- 

In the above we thus found the general solutions of the Schrodinger equation for the time-dependent linear potential, 
which corresponds only to the quadratic-form invariant (3.2). It is concluded here that the solutions obtained above 
does not form a complete set of solutions of this time-dependent Schrodinger equation, and that Bekkar et al.'s solution 
and my solution presented here will constitute together such complete set of solutions of the Schrodinger equation. 

IV. DISCUSSIONS AND CONCLUSIONS 

(i) In the present paper we show that since Bekkar et al.'s solution [3] has not yet contain those corresponding to 
the quadratic invariant, it is not the true general solution of the Schrodinger equation for the time-dependent linear 
potential. Instead, it is the solution corresponding only to the linear L-R invariant. The obtained solution here is the 
one that corresponds to the invariant (3.2), which is of the quadratic form. Since the linear and quadratic invariants 



\n,t) 



Schr 



exp 



J0„(t) 
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form a complete set of L-R invariants, Bckkar et aVs solution [3] and my solution presented here constitute such 
complete set of solutions of the Schrodingcr equation involving a time-dependent linear potential. 

(ii) It is well known that in quantum optics there are three kinds of photonic quantum states, i.e., Fock state, co- 
herent state and squeezed state. From my point of view, the calculation of the variations of creation and annihilation 
operators (aJ, a) of photons under the translation (e.g., V\ of (3.4)) and squeezing transformation (e.g., V% of (3.4)) 
operators shows that the variations of a' and a are exactly analogous to that of space-time coordinate variations un- 
der the translation, Lorentz rotation (boosts) and dilatation (scale) transformation [9] and thus these three quantum 
states (coherent, squeezed and Fock states) correspond to the above three conformal transformations, respectively. I 
think that this connection between them is of physical interest and deserves further consideration. 

(iii) Guedes recently stated that in order to obtain the general solutions of the Schrodinger equation one must follow 
the L-R invariant theory step by step [8]. I don't approve of this point of view, however. Personally speaking, in fact, 
the L-R method has only one step, namely, the particular solution of the time-dependent Schrodingcr equation is 
different from the eigenfunction of the invariant only by a time- dependent phase factor. In Ref. [3] and [5], although 
we follow the L-R method step by step, what we obtained still cannot be viewed as the general solutions of the 
Schrodinger equation. For this reason, I think that "step by step" is not the essence of getting the general solutions 
of Schrodingcr equation. Instead, the key point for the present subject is that one should first find the complete set of 
all L-R invariants of the time-dependent quantum systems under consideration. For some systems in the Hamiltonian 
there may exist no such closed Lie algebra as (1.1), the complete set of exact solutions can be found by working in 
a sub-Hilbert-space corresponding to a particular eigenvalue of one of the invariants, namely, only in the sub-algebra 
(quasi-algebra) corresponding to a particular eigenvalue of this invariant will such time-dependent quantum systems 
(which have no closed Lie algebra) be solvable [10]. For the timc-dependent quantum systems, there are no other 
eigenvalue equations of Hamiltonian than that of the L-R invariants with time-dependent eigenvalues. The complete 
set of invariants, instead of the time-dependent Hamiltonian, can describe completely the time-dependent quantum 
systems. For this reason, it is essential to find the complete set of invariants for the time-dependent Hamiltonian of 
a given quantum system. 

(iv) In the Ref. [2] , the author says that to the best of his knowledge there was no publication reporting the solution 

2 

of the Schrodingcr equation for the system described by H(t) = ^ + f(t)q without considering approximate and/or 
numerical calculations. I think that this is, however, not the true case. In the literature, at least in the early of 
1990's, Gao et al. had reported their investigation of the driven generalized time-dependent harmonic oscillator which 
is described by the following Hamiltonian H(t) = \[X(t)q 2 + Y(t)(pq + qp) + Z(t)p 2 } + F(t)q [4]. It is believed that 
my solution presented here is only the special case of what they obtained [4]. 
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